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We demonstrate that the clustering statistics and the corresponding phase transition to non- 
equilibrium clustering found in many experiments and simulation studies with self-propelled particles 
(SPPs) with alignment can be obtained from a simple kinetic model. The key elements of this 
approach are the scaling of the cluster cross-section with the cluster mass - characterized by an 
exponent a - and the scaling of the cluster perimeter with the cluster mass - described by an 
exponent f3. The analysis of the kinetic approach reveals that the SPPs exhibit two phases: i) an 
individual phase, where the cluster size distribution (CSD) is dominated by an exponential tail that 
defines a characteristic cluster size, and ii) a collective phase characterized by the presence of non- 
monotonic CSD with a local maximum at large cluster sizes. At the transition between these two 
phases the CSD is well described by a power-law with a critical exponent 7, which is a function of a 
and j3 only. The critical exponent is found to be in the range 0.8 < 7 < 1.5 in line with observations 
in experiments and simulations. 

PACS numbers: 87.18.Gh, 05.65.+b, 87.18.Hf 
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I. INTRODUCTION 

Many experimental self-propelled particle (SPP) sys- 
tems JlU3| , from actin filaments driven by molecular mo- 
tors U1J gliding and swimming bacteria > to active 
colloidal particles @, 9] exhibit a remarkably rich clus- 
ter dynamics. In particular, the formation of large mov- 
ing polar clusters has been observed in most of these 
examples. Such clusters are formed by particles that 
move roughly in the same direction. Simulations of 
simple models of SPPs have revealed that similar clus- 
tering dynamics are observed for SPPs with both, po- 
lar [ll]l - fl2j and nematic fl3l - [T5| alignment interactions. 
More detailed models display similar clustering statis- 
tics, e.g., in simulations of self-propelled rods 
self-propelled disks [Til , particles with hydrodynamic 
coupling [2l|, and swimming particles with flagella, e.g., 
in sperm cells (22j . 

In many of these systems, a peculiar phase, character- 
ized by the existence of remarkably large moving clusters, 
has been observed. One refers to this phase of collec- 
tive motion as non-equilibrium clustering. This phase 
appears often as an intermediate phase between a com- 
pletely disordered phase, with homogeneous density, and 
a phase with global orientational order. An experimen- 
tal case (myxobacteria) and a simulation example (self- 
propelled rods) are shown in Fig. [TJ The occurrence of 
large clusters is strongly correlated with a crossover in 
the shape of the cluster size distributions (see Fig. \Dp, 
d). Figure [2] illustrates recent findings in simulations of 
self-propelled hard rods [H, HH an d in simulations with 
a modified Vicsek model with nematic alignment [l3l - [T5| . 
The collective clustering phase is characterized by a non- 
monotonic cluster size distribution (CSD) with a charac- 
teristic peak at large cluster sizes, see CSD for high den- 




FIG. 1: Examples of non-equilibrium cluster formation in self- 
propelled particle systems, (a) Large moving clusters in ex- 
periments with myxobacteria, where it has been found that 
the cluster size distribution (CSD) is a function of the cell 
density as shown in (b) (see [6|]). (c) Snapshot of moving 
clusters in self-propelled rod simulations, where the CSD is 
known to be function of the particle density, as illustrated in 
(d) (see [3). 



sities in Fig. [T] At the onset of the collective clustering 
phase, the CSD follows a power law with a characteristic 
exponent. It has been observed in SPP experiments 0,0] 
and simulations 0, [ItJ that the CSD can be power-law 
distributed with the exponent typically varying between 
0.85 and 1.35. While the range of exponents indicates the 
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FIG. 2: Examples for the onset of the collective (clustering) 
phase (blue) and global order (orange) in simulations, (a) 
Self-propelled hard rods in two-dimensions with asp ect ratio 
k = 13, v denotes the fraction of covered area [231. 1241]. (b) 
Modified Vicsek model with nematic alignment with density 
p = 0.25 [lij. The noise value for the transition to global 
order was taken from simulation with > 10 6 self-propelled 
particles [3] , while the noise value for the onset of clustering 
in a system with 16.384 self-propelled particles is based on 
the results of [l5| . 



absence of a universal scaling at the onset of the collec- 
tive clustering phase, it is not yet clear what determines 
the actual value of the exponent. 

The examples displayed in Figs. [1] and [2] indicate that 
the collective clustering phase often occurs in the ab- 
sence of global order. Other models, e.g. the classical 
Vicsek model displa y a similar phase, but in the pres- 
ence of global order [ill . H2I l25j . A word of caution is in 
place here. Most investigations that focused on the clus- 
ter statistics (both, in experiments and simulations) were 
carried out with small to intermediate numbers (around 
10 2 — 10 3 ) of self-propelled particles, while investigations 
regarding the onset of global order were conducted with 
rather large numbers (> 10 5 ). Hence, it is not completely 
clear if and how the system size, respectively the parti- 
cle number, affects the onset of this collective clustering 
phase. 

Here, we focus on the emergence of the collective clus- 
tering phase in the absence of global order, i.e. we assume 
that cluster velocities are uncorrelated. We propose a ki- 
netic clustering theory, based on the approach presented 
by us in [HI ] , and extend our earlier treatment to coagula- 
tion and fragmentation kernels that depend, respectively, 
on the scaling of the cluster cross-section with cluster 
"mass", i.e. number of particles in the cluster, and the 
scaling of cluster perimeter with cluster mass. The scal- 
ing of the cluster cross-section and cluster perimeter are 
characterized by exponents a and /3, respectively. Fur- 
thermore, we present a comprehensive analysis of this 
kinetic cluster model, including a finite size (FS) study 
of it. The FS study of the kinetic model reveals that: i) 
the transition to the collective clustering phase - charac- 
terized by a non monotonic CSD - is generic to all SPPs 
with either a polar or nematic alignment mechanism, and 
ii) that there exist a critical asymptotic CSD exponent 



7 which is function of a and (3 only. By studying the 
physically meaningful parameter space a-/3 we find that 
7 always falls in the range 0.8 < 7 < 1.5. Notice that in 
the absence of cluster-cluster correlations, this transition 
marks the onset of collective motion as observed in [6j. 
Though the simplified clustering theory is strictly speak- 
ing only valid in the absence of global order, we show 
through simulations that a comparable picture holds for 
systems where global order is observed. 

The paper is organized as follows. In section [H] we 
present our simple kinetic clustering model under the as- 
sumption that there are no cluster-cluster correlations 
and show that for any finite system there are two clus- 
tering phases. We perform a system size analysis of the 
kinetic model in Sec. IIII1 where we report on the scaling 
properties of the obtained cluster size distribution (CSD) 
which we find to depend on the scaling of the cluster 
cross-section and cluster perimeter. In Sec. IIIII we study 
cluster formation in a SPP model that is known to exhibit 
cluster-cluster correlations and contrast the obtained re- 
sults with our cluster-cluster uncorrelated clustering the- 
ory The implications and limitations of the proposed ki- 
netic approach are discussed in Sec. El 



II. A KINETIC MODEL FOR CLUSTERING 

We look for a description of the clustering process in 
terms of (ni(i)), (ri2(t)}, (nAr(i)), where (n,i(t)) repre- 
sents the average value at time t of clusters formed by i 
particles, and N the number of particles in the system. 
To ease the notation, we refer to (rii(t)) as Ui{t). Our 
strategy consists on deriving a coagulation Smoluchowski 
equation with fragmentation for these rii(t) objects. For 
an introduction to this kind of coagulation equations we 
refer the reader to [26]. In SPP systems, the only con- 
served quantity is the number of particles N. Neither 
the overall orientational order nor the number of clus- 
ters are conserved. Thus, our generalized coagulation 
Smoluchowski equation should conserve the number of 
particles. We start by simplifying the clustering dynam- 
ics. Since we are dealing with self-propelled particles, we 
assume that clusters are the result of either an explicit 
or an effective alignment mechanism such that particles 
inside a cluster move coherently in the same direction. 
This implies that clusters move at speeds comparable to 
that of the individual particles, independently of the size 
of the cluster. Below we also discuss what can be ex- 
pected if this condition is relaxed. On the contrary, clus- 
ters of passive particles, driven by thermal fluctuations, 
are such that their mobility decreases as function of their 
size. We stress that the assumption of a size-independent 
cluster speed evidences the non-equilibrium nature of our 
simple cluster theory. We further assume that (binary) 
collisions among clusters may result in cluster-cluster fu- 
sion, and neglect the possibility of cluster fragmentation 
induced by cluster-cluster collision. This assumption is 
justified as long as the dynamics is overdamped, and is 
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in accordance to what is observed in experiments with 
gliding bacteria [6], and simulations with self-propelled 
rods [16]. We simplify the cluster fragmentation dy- 
namics by assuming an evaporation-like process by which 
clusters shrink in size by loosing one by one those par- 
ticles that are on the cluster boundary. More complex 
fragmentation process can be ignored as exponentially 
unlikely events. In summary, we are assuming an irre- 
versible clustering process in which a cluster of mass j 
can undergo the following "reactions" : 



Cj + Ck Cj+fc, Cj -4 Cj-i + C\ 



(1) 



Notice that since the reaction Cj+fc —> Cj + Ck does not 
occur, the process is, in this sense, irreversible. Now, 
we look for a description of the process in terms of the 
(average) number rij(t) of clusters with j particles (i.e., 
number of Cj's) at time t, whose time evolution takes the 
form: 
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where the dot denotes the time derivative, Bj represents 
the rate at which a cluster of mass j looses particles, 
defined as 



^ ~ cP 3 



(3) 



and Aj k is the collision rate between clusters of mass j 
and k, defined by 



l 3,k 



L 2 



(4) 



where L 2 is the area of the two-dimensional space where 
particles move and v c is the cluster speed, which we as- 
sume to be v c ~ vq with vq the speed of individual par- 
ticles. Notice that Eqs. ([2]) are such that ^mn m = 0, 
and thus the number of particles N = ^ m n m (t) is con- 
served. In Eq. ([3]), d 2 /D is the typical time a particle 
located on the cluster boundary needs to detach from a 
cluster, with d the maximum distance two particles can 
be apart to be still considered as connected and D the 
diffusion coefficient with respect the center of mass of 
the cluster. We stress that in Eq. @ we have assumed 
that cluster are uncorrelated. Several statistical features 
of the model can be given explicitly, if one assumes that 
SPP follow a dynamics where their translation is deter- 
mined by a constant speed and the direction of motion is 



subject to alignment interactions with neighboring SPP 
and an angular noise with amplitude r\. An example of 
such dynamics is given by the following equations of mo- 
tion: 



9t+i 

j 



arg 



(5) 
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where j is the index of the particle, 6j defines (in two 
dimension) the particle moving direction, x*- denotes the 
position at time t, the sum is taken over all particles 
within a unit distance of j, and 77* is a uniformly dis- 
tributed random variable such that — rj/2 < 77* < t]/2, 

(77*-) = 0, and (rfjrfl) = {rj 2 /12)S t ,t'Sj, k - The dynamics 
given by Eq. (fp is frequently used in simulation stud- 
ies, e.g. in [ToPl^l . l27l. I2 H (for a review, see [l|). More 

specifically, f(e l6k ,e l9 ^ — e l6k defines a polar align- 
ment rule as in the Vicsek model 
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sign[cos(#). — 0j)]e k defines SPP with a nematic align- 
ment rule as used [H, [13] ■ Assuming that the SPP obey 
Eq. ||SJ), the diffusion around the center of mass of a polar 
cluster of SPP can be approximated, as detailed in (l5j . 
by D = i7g[l-(sin(ri/2)/(7i/2)) 2 ]A£, where v is the speed 
of individual particles, and At sets the discrete time step 
for a typical SPP dynamics such as the one given by 
Eq. ([5]). By making a Taylor expansion we find that 
D ~ VQ/l2i] 2 At. Notice that for non-interacting random 
walkers subject for the same kind of noise, D cx v^/r] 2 . 
In addition, we have to consider that the splitting rate 
Bj has to be proportional to the number of particles 
on the cluster boundary. In Eq. ([3]), we have assumed 
that the perimeter I of a cluster scales with its mass j 
as I oc jP, with /3 constrained to 1/2 < ft < 1. Let 
us provide a physical context to this assumption, always 
assuming we are in a two dimensional space. Surface 
tension would tend to minimize cluster perimeters, and 
so clusters would be round and (3 = 1/2, as observed in 
liquid- vapor drops 29]. On the other hand, for random 
deposition of particle as in classical percolation, w 1. 
If clusters are "chains" of particles then also /3 = 1. In 
short, we expect 0.5 < /? < 1. 

The collision rate Ajk is derived in analogy to the col- 
lision rate in kinetic gas theory [30] , which means we 
assume clusters move ballistically in between collisions, 
which leads to a collision rate proportional to the sum 
of the scattering cross section of the involved clusters, 
the relative cluster speed, and cluster density. Due to 
the latter assumption, Aj t k oc L~ 2 . The relative cluster 
speed can be assumed to be not too different from that of 
the individual clusters, v c . In SPP systems, we find often 
that individual cluster speed is close to that of individ- 
ual particles, i.e., v c ~ vq. If v c exhibits a dependency on 
the cluster size of the form v c oc the exponent £ can 
be absorbed into a. For simplicity, here we assume that 
v c ~ vq. For instance, the speed of SPP clusters can be 




FIG. 5: Exponents 7 as function of the exponents a and /3. 
In (a) a — /3, while in (b) a = 0.5 and /3 > a. The inset in 
(b) shows that when /3 > a, £ < 1; for a = /3, £ — 1, see Eq. 
CEO]). 



FIG. 3: For a given system size N, the cluster size distribution 
(CSD) p(m) is monotonically decreasing for P < P c (a), while 
for P > P c the distribution exhibits a peack at large cluster 
sizes (b). At P = P c the distribution p(m) is a power-law. 
The transition between these two behaviors, (a) and (b), can 
be represented in a more phase transition-like form, by looking 
at first two moments of the distribution p(m), (c) and (d), 
respectively. Other parameters: N — 512 and f3 = 1/2. 
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FIG. 4: (a) Scaling of p(m) at the critical point P C (N) for 
various system sizes N. Notice that for N — > 00 the distribu- 
tion is given by Eq. ((8}. (b) Below and above P C (N), p(m) 
exhibits the scaling given by Eq. (|SJ). (a) and (b) corresponds 
to a = P = 1/2. For the value of 7 see Fig. [5] 
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approximated by v c = vosm(rj/2)/(r]/2) oc vo{l — r/ 2 ) [l5j . 
The scattering cross-section of a 2D cluster is the (av- 
erage) projection of the cluster on a given axis. Thus, 
scattering cross section has to be always less or equal 
to the cluster perimeter. If the scattering cross section 
of a cluster scales as oc (Tq j a , the previous observation 
implies that a < ft is the only meaningful physical sce- 
nario. Below, we will see that there are two qualitatively 
very different physical scenarios, a = j3 and a < /3. Fi- 
nally, q s represents the probability that a cluster-cluster 
collision resulted in a successful fusion of clusters. Here, 
we assume it is a constant but certainly it can depend on 
various (intensive) variables of the actual system. In par- 
ticular, it depends on the alignment symmetry. A very 
rough assumption would be to assume that if q s = qo for 
the ferromagnetic alignment, q s = qo/2 for the nematic 
one. Eqs. @ are scaled and transformed into a dimen- 
sionlcss form by dividing Eq. <j3j) and (J4j) by D/d 2 , which 
leads to the following dimensionless parameter: 



P 



q s a a d 2 v c 
L 2 D 



(6) 



Thus, Eq. © reduces to Bj = f , and Eq. @ to 
Aj.k = P(j a + k a ). The parameter P controls the rela- 
tive weight of fragmentation with respect to coagulation. 
For model type considered here, we obtain P oc (Lrj)~ 2 . 
The computation of this critical point P can be done by 
directly studying the stability of the individual phase, 
and ignorin g th e actual shape of the CSD, as recently 
proposed in [31j . 

From direct numerical integration of Eqs.©, in its di- 
mensionless version, with initial condition rij(t — 0) = 
NSj.i, we find that the (weighted) cluster size distribu- 
tion (CSD), defined as 

p(m,t) 



mn m (t) 
N 



(7) 



reaches a steady state, i.e., p(m,t — >• 00) = p(m). This 
also implies that the number of cluster in the systems, 
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M (t) = NJ2 m ~ 1 Pm{t), also reaches a steady state value 
as t — > oo. In the literature the term CSD is frequently 
used to refer to p(m, t) — n m (t)/M(t). The disadvantage 
of this definition is that its normalization constant, M(t), 
varies with time. While p(m) refers to the probability of 
finding a cluster of size m, p(m) indicates the probabil- 
ity of a randomly selected particle to be in a cluster of 
size m. Given a and /3, the CSD p(m) depends on N 
and the value of P. Fig. [3] summarizes the clustering 
behavior with P for a given N. Notice there exists a 
critical P C (N) that separates two different clustering be- 
haviors. For small values of P, the distribution p(m) is 
monotonically decreasing, dominated by an exponential 
tail. As P — > P c , p(m) approaches a power-law, with a 
system size cut-off. For P > P c , p(m) is non-monotonic 
and exhibits a peak at large cluster sizes emerges. This 
dramatic change of behavior at P c unveils a phase tran- 
sition, which is evidenced by Fig. [21(c) and (d) that show 
( m ) = ^rnp(rn) and a 2 — XX m — (m)) 2 p{m), respec- 
tively. Below, we will see that this transition is a genuine 
phase transition - in the thermodynamical sense - only 
for a — (3. 

III. SCALING PROPERTIES 

For every N, we numerically estimate P C (N) as the 
point at which p(m) is no longer monotonically decreas- 
ing. At the critical point p c (N), p(m) scales as: 

p(m; P c ) oc m -7 , (8) 

where 7 is a critical exponent that depends on a and 
(3. This fact is illustrated by Fig. |3{a) that shows that 
pirn; P c ) follows the scaling given by Eq. (|SJ up to a given 
cluster size m*(N) above which the finite size of the sys- 
tem becomes evident. We find the in general, below the 
critical point, the following scaling is obeyed by p(m): 

p(m; N, P, a, 0) = K p(m; bN, b -e P, a, (3) , (9) 

where K is a constant, see collapse of the curve in 
Fig. Hb). For a = (3, K = 1, while for a ^ /3, K 
depends on b, N, and j3. This means that if we double 
the system size and reduce P by 2 - ^, we fall on the same 
distribution p(m) for a = /3, and for a < j3 there is a 
multiplicative constant K. For P < P c , we find that the 
CSD is well fitted by p(m) ~ m -7 exp(— m/rh(P)). 

For P > P c , we observe that there is a shift with N 
of the peak that emerges at large cluster sizes, Fig. [3Jb) 
The scaling given by Eq. © works up to a given size 
m m (N) which is given by the minimum of p(m). We are 
interested in knowing the behavior of the peak with the 
system size N. In order to answer these questions we 
study the scaling of Q(N) = Ylm N p( m )- We find that 
for a = (3, Q(N — > 00) — > 1. This indicates the presence 
of a collective phase in which most of the particles are 
part of large clusters. Notice that this does not mean the 
formation of a single giant cluster; the width of the peak 



does not shrink to zero. For a < (3, Q does not converge 
to 1 as N — > 00. This is because for a < (3, in the limit 
of N — > 00 the transition does not occur. Nevertheless, 
for any finite N, we always observe the transition to the 
collective phase. This observation can be understood by 
studying the system size behavior of p c (N). We find that: 

P C (N) oc N~t , (10) 

with £ = 1 for a = /3, and £ < 1 for a < (3. For a fixed 
value of a, £ becomes a function of f3 as illustrated in the 
inset of Fig. [5] (b) for a = 1/2. In the thermodynam- 
ical limit, with N —> 00 and L —> 00, while p = N/L 2 
constant, we can estimate the critical point. Combining 
Eqs. © and ([TO]) we obtain, for a given set of parameters 
D, Vq, (To, and d, the critical density p c above which a 
peak at large cluster sizes emerges: 



From Eq. (jlll) is clear that when £ < 1, p c diverges with 
the system size N. Only for £ = 1, p c is a finite quantity 
in the thermodynamical limit. Thus, for a = (3 there 
is a critical density p c above which the collective tran- 
sition occur. On the contrary, for a < j3, p c — > 00 and 
any finite densitiy po corresponds to the mono-disperse 
phase. Eq. (jTTJ) predicts the critical noise intensity r) c 
below which one expect the collective clustering phase 
to occur as r\ c oc \J pN^^ 1 , from which follows that for 
£ = 1, the critical noise scales as i] c oc ^fp. Interestingly, 
a similar scaling was reported for the onset of collective 
motion in simulations with the Vicsek model [l(J HH HI] . 

IV. CLUSTERING IN THE PRESENCE OF 
GLOBAL ORIENTATIONAL ORDER 

The simple set of Smoluchowski-type Eqs. ([5} has 
been derived assuming that clusters move in an uncor- 
related fashion, as observed for instance, in bacterial 
experiments @. However, several SPP models exhibit 
long-range orientational order and macroscopic struc- 
tures such as bands, which indicates that there are strong 
correlations among cluster velocities. Thus, a priori 
we cannot expect Eq. © to describe quantitatively the 
cluster dynamics of such systems deep in phases with 
strong cluster-cluster correlations. Interestingly, it has 
been already observed that the Vicsek model (VM) in its 
(global) ordered phase displays CSDs that can be power- 
law distributed, with 0.8 < 7 < 1.3, according to [1 ll |25(| . 
Here, we carry out simulation in the Vicsek model (VM), 
i. e. Eqs. ([5]) with a choice of / = e l6k . Our simula- 
tion results confirm that indeed the VM exhibits power- 
law distributed CSDs with exponents in this range, and 
show that this occurs close to the well-known disorder- 
order transition of the VM. Furthermore, we find that 
the cluster statistics in the VM, below and close to the 
critical point - i.e., when cluster-cluster correlations are 
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FIG. 6: Simulations snapshots of the VM for various noise intensity values. Particles belonging to the same cluster are labeled 
with the same color. As it can be seen in Fig. [7] for r\ — 1.0 the system is disordered, while for r\ < 1 there is preferred direction 
of motion. Though the snapshots for 77 = 1.0 and 77 — 0.9 look comparable, their CSD is remarkably different: an exponential 
distribution for the former, and a power-law for the later. Traveling bands in the VM are composed of many correlated clusters 
(77 = 0.7). Finally, for very low 77 values most particles form part of very large clusters (77 = 0.05). 



(a) (b) 




FIG. 7: Clustering in the VM. (a) shows the CSD for various 77 
values, while (b) displays the scaling close to the critical point 
?7 C . The average (normalized) cluster size (m)/N as function 
of 77 is shown in (c), and the orientational order parameter 
4> vs. 77 in (d). While the system is disordered, the CSD is 
dominated by an exponential tail. Below the critical point, 
the CSD is a power-law with an exponent 7(77). Close to the 
critical point 7 = 1.3 and as 77 is decreased 7 approaches 1. 
For very low 77 values, the CSD decreases even slower than 
1/m (see Fig. [6}. 



weak or absent - resembles that obtained with the kinetic 
clustering model. As expected, far away from the critical 
point and well in the ordered phase, cluster-cluster cor- 
relations induce effects that cannot be accounted by the 
kinetic clustering model. In short, the simulation data 
strongly suggests that close to the disorder-oder critical 
point, the VM exhibits a transition to a collective clus- 
tering phase as the one described with the kinetic model. 

The order-disorder transition in the VM is observed 
when the noise intensity 77 is decreased below a crit- 
ical value 77*. Orientational order is characterized by 
4> = I J^j ex P(^fc)l/^- I n the disorder phase, the clus- 



ters cover homogeneously the space (Fig.[5J) and the CSD 
is dominated by an exponential tail, see Fig. [7J If 77 is 
decreased below 77*, velocity- velocity correlations among 
particles become important, <j) increases, and clusters 
grow significantly in size. Below the critical point p*, 
the CSD is a power-law with an exponent in the range 
[0.8,1.3], which results also in a change of behavior of 
(m)/N , as shown in Fig. [7J The power-law distributed 
CSD indicates that the system displays arbitrary large 
clusters. Indirectly, this also means that particle-particle 
correlations are long-ranged, i.e., arbitrarily large as the 
cluster sizes. Notice that this does not necessary imply 
the existence of cluster-cluster correlations. Nevertheless, 
in the VM at low values of the noise intensity cluster- 
cluster correlations become evident. Arguably, due to 
these cluster-cluster correlations, we observe a clear devi- 
ation from what the theory predicts for 7/ « 77* . More- 
over, in the VM, 7 is a function of 77. We observe that 
close to 77*, 7 ss 1.3, and as 77 is decrease, 7 approaches 
1. It is close to 77* that we observe the emergence of a 
traveling band. Notice that this band is not a connected 
component, but a cloud of highly correlated moving clus- 
ters. We find that the CSD is power-law distributed in 
the band regime. For 77 << 77* the system is highly or- 
dered and clusters move roughly in the same direction, 
which means that cluster collisions are less frequent. On 
the other hand, low values of 77 also imply that spreading 
of a cluster around its center of mass is also small. In 
summary, at low noise values the cluster dynamics slows 
downs significantly in the VM, which results in a non 
monotonous response of (m)/N with 77, Fig. [7J Interest- 
ingly, we observe that at extremely low 77 values the CSD 
decreases even slower than 1/m and (m)/N experiences 
a sharp increase. 
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V. CONCLUDING REMARKS 

We have argued that SPP systems exhibit for a finite 
system size two phases: a mono-disperse and a collec- 
tive clustering phase, characterized by a non monotonic 
CSD. Assuming that the moving directions of clusters are 
uncorrelated, Eq. justifies the existence of these two 
phases. At the transition point, the CSD is a power-law 
characterized by an exponent 7 that depends, according 
the proposed kinetic clustering model, on the scaling with 
cluster size of the cluster cross-section - characterized by 
an exponent a - and cluster perimeter - characterized by 
an exponent j3. A systematic study of the a— f3 parameter 
space revealed that 7 always falls in the range between 
0.8 and 1.5, which is consistent with experimental obser- 
vations [H, Q and simulations 0,113] • However, we learn 
that in the thermodynamical limit, only for a = f3 the 
above mentioned transition occurs, while otherwise, the 
system remains in the mono-disperse phase. It is worth 
pointing out that for the special case a — (3 the kinetic 
model predicts a critical density value that is indepen- 
dent of the system size. Hence, the properties for the 
transition to nonequilibrium clustering found in the ki- 
netic model for intermediate values of iV hold true in the 
thermodynamic limit. 

In addition to the study of the kinetic clustering theory 
given by Eq. ((2J, we characterized the cluster statistics 
of the Vicsek model [l(| • Despite the fact that cluster- 
cluster correlations and global order emerge in the sim- 
ulations of the Vicsek model that are not treated in 
Eq. ([2]), we found that a transition to a collective cluster- 
ing phase is also present in the VM. Moreover, 7 falls in 
the expected range if simulations are performed with fer- 
romagnetic alignment, as shown here and in [ill l25|. as 



well as with nematic alignment [13 . Notably, we found 
that the CSD is power-law distributed for a range of 
noise intensity values r\ close to the order-disorder critical 
point, with 7 function of r\. In summary, the existence of 
two clustering phases is found as well in SPP systems in 
the presence or absence of global orientational order. 

Our kinetic approach is complimentary to the ongoing 
effort of deriving coarse-grained nonlinear field equations 
for the description of active matter |32h39| | (for a review, 
see @, EO]) that are better suited to capture large-scale 
structures or the emergence of global orientational or- 
der in large systems. In contrast, the observed clustering 
phenomena are typically studied in the absence of long- 
range order, on a small scale, and for intermediate system 
sizes. It is also worth noting, that the non-equilibrium 
clustering phase leads to (apparent) giant number fluc- 
tuations [6| which often have similar properties as the 
ones predicted by Toner and Tu |32j in the phase with 
long-range orientational order. 

We expect a similar transition and cluster dynamics as 
the one described here for all SPP systems where particle 
speed is not strongly affected by the local density. The 
presence of a density dependent speed can dramatically 
change the above given clustering picture, since now large 
size clusters are prone to slow down significantly. As 
result of this effect, the SPPs can form a single large 
cluster which coexists with a background gas of particles 
in the individual phase, as found in [4ll - [44| . The cluster 
dynamics of these systems is likely to be related to an 
equilibrium-like phase separation, as suggested in [4ljj . 
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